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ABSTRACT 



^ Aims. The purpose of this work was to obtain diffusion coefficient for the magnetic angular momentum transport and material transport in a 
rotating solar model. 

Methods. We assumed that the transport of both angular momentum and chemical elements caused by magnetic fields could be treated as a 
diffusion process. 

CO Results. The diffusion coefficient depends on the stellar radius, angular velocity, and the configuration of magnetic fields. By using of this 
(^-^ , coefficient, it is found that our model becomes more consistent with the helioseismic results of total angular momentum, angular momentum 
■ density, and the rotation rate in a radiative region than the one without magnetic fields. Not only can the magnetic fields redistribute angular 
QO momentum efficiently, but they can also strengthen the coupling between the radiative and convective zones. As a result, the sharp gradient of 
the rotation rate is reduced at the bottom of the convective zone. The thickness of the layer of sharp radial change in the rotation rate is about 
0.036 Rq in our model. Furthermore, the difference of the sound-speed square between the seismic Sun and the model is improved by mixing 
the material that is associated with angular momentum transport. 
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1. Introduction 

In standard stellar models, it is assumed that there are no rota- 
tion and magnetic fields. Although this framework is strongly 
supported by helioseismology, recent advances in the study of 
solar structure shows that differences exist between the Sun 
and this model. Those differences are small and display a very 
interesting systematic behavior, which is far from satisfactory 
(Christensen-Dalsgaard |20()2] l. The most striking difference is 
the bump in the sound speed just beneath the convection zone. 
Noting both cc T/fi and the //-gradient caused by micro- 
scopic diffusion and element settling in the radiative region, 
Gough et al. (119961 ) and Christensen-Dalsgaard et al. ( 119961 
12002 1) argued that material mixing may be important but has 
been neglected. 



Helioseimology also shows that rotation is almost uniform 
in the solar radiative region: the Sun has a slow rotation core, 
the angular velocity Q has a small radial gradient and large 
latitudinal gradient in the convective zone, and the radial an- 
gular velocity gradient is positive at low latitudes and negative 
at high latitudes in the tachocline (Schou et al. 1998; Chaplin 
et al. 119991 ). But the calculations of rotational models show 
that the sun has a fast rotation core and a large gradient of the 
rotation rate in the radiative region (Pinsonneault et al. 119891 
Charboyer et al. 119951 ). which disagrees with the helioseismic 
results. Thompson et al. (2003') argues that some mechanisms 
of angular momentum transport may be missed in these mod- 
els. Magnetic and gravity waves mechanisms have been pro- 
posed; however, what mechanisms are causing the nearly uni- 
form rotation in the radiative region are poorly understood. 
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Rotation is a property that virtually all stars possess. 
Rotation affects stellar structure mainly in two ways. The first 
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is an immediately dynamical effect caused by centrifugal ac- 
celeration of the hydrostatic balance. In most cases, however, 
this effect is modest and can be included through modifying 
the equations of stellar structure. Although the equations of 
stellar structure of a rotating star are three-dimensional, the 
method of Kippenhahn & Thomas (1 19701 1. which uses the mass 
contained within an equipotential surface, M^, as an indepen- 
dent variable, allows a one-dimensional evolution code to be 
modified to incorporate the hydrostatic effects of rotation. This 
method was modified by Endal & Sofia ( 1976) and Meynet & 
Maeder ( |1997l l to apply to shellular rotation (Zahn |1992| l, a 
rotation rate that depends only on the radial coordinate r to 
a first approximation. The other way is an effect arising from 
the redistribution of chemical elements due to the instabilities 
caused by rotation. This effect is much more important and it 
has been studied by many investigators (Endal & Sofia 1978. 
Pinsonneault et al. 1989; Chaboyer & Zahn [19921 Zahn [19921 
Meynet & Maeder tl997j Maeder 1997; Maeder & Zahn[T998l 
Maeder & Meynet HOOOl Huang 2QP4a 2004b ). The method 
of Kippenhahn & Thomas ( 1 19701 ) as modified by Meynet & 
Maeder ( 1997 ) is used in our model. 

Magnetic fields are another property of stars. Their effects 
on a star create a complicated problem. The magnetic field is 
involved in most problems in astrophysics. But the generation 
of magnetic fields is still a controversial problem. They may 
be fossil fields (CowHng [l945l Moss [T987i Braithwaite 2004), 
which are remnants of the star's formation, or may be generated 
by a convective stellar dynamo (Parker [19791 Charbonneau & 
MacGregor 12001b in the convective zone or generated by a 
Tayler-Spruit dynamo (Pitts & Tayler ^1986', Spuit 2002) in a 
differential rotation star In this paper we assumed that the mag- 
netic fields exist. We only consider the process of magnetic 
angular momentum transport and material mixing. Although 
the magnetic angular momentum transport has been treated by 
many investigators (Charbonneau & MacGregor 119921 119931 
SpruitlSOOS! Maeder & Mevnet i2003ll2004b . it is still an open 
question. 

In this paper, we focus mainly on magnetic angular mo- 
mentum transport and the mixing of elements associated with 
angular momentum redistribution. In Sect. 2 we give the equa- 
tions of stellar structure of a shellular rotation star. In Sect. 3 
we deduce the diffusion coefficient for magnetic angular mo- 
mentum transport and material mixing that is due to angular 
momentum redistribution. Then, in Sect. 4 we give the results 
of the numerical calculation. We then discuss our results and 
conclude in Sect. 5. 



isobar surface is denoted by S p, and the volume enclosed by 
the isobar surface by Vp. For any quantity q, which is not con- 
stant over an isobar surface, a mean value is defined by 



< q >- 



r! 



qdcr , 
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where dcr is an element of the isobar surface. 

The equations of stellar structure of a shellular rotation star, 
which were developed by Kippenhahn & Thomas d 19701 1 and 
Meynet & Maeder (,1997 ), are as follows: 
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where rp is the radius of a sphere enclosing the volume Vp, i.e., 
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The Mp is the mass inside the isobar, and 
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Here ge is the effective gravity. The a in Eq. © is a scalar 
Equation Q is the hydrostatic equilibrium equation including 
centrifugal force, while Eq. ^ is the equation of energy trans- 
port under the effects of rotation. The nondimensional rotating 
corrective factors fp and fj depend on the shape of the isobars. 

Assuming that the shapes of isobars are spheroids with 
semi-major axis a and semi-minor axis b, and given the angular 
velocity distribution, using definitions ([7]i and ([TJ, i.e.. 



2. Equations of stellar structure of a shellular 
rotation star 

The stars whose angular velocity is constant on their isobars 
are called shellular rotation stars (Meynet & Maeder 1997 ); the 
isobar surfaces are given by (Meynet & Maeder |1 9971 1 



>!',, = (J) -H ]-Q}p-sirp-6 - constant. 



(1) 



where <1) is the minus gravitational potential, Q. the angular ve- 
locity, r the radius, and 6 the colatitude. The area of such an 
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we can get the values of a and b for any given Mp and r^; i.e., 
the surface of the isobar is determined. Thus, the average values 
of the effective gravity and its inverse (< g^ > and < g^^ >) 
can be obtained. Then the fp and fj are also obtained. In the 
computation, the Roche mode (Kippenhahn & Thomas [19901 1 
was used to compute the gravitational potential. 
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3. Angular momentum transport and the mixing of 
elements 

3.1. Hydrodynamical instabilities 

The mechanisms that redistribute angular momentum and the 
chemical elements in a rotating star can be divided into two 
categories according to the time scale involved. 

The first category is that of the dynamical instabilities, oc- 
curring on the dynamical timescale. If the dynamical unsta- 
ble gradient occurs in star, it can be instantaneously smoothed 
(Endal & Sofia |1978] ). One of the dynamical instabilities is the 
convective instability. We suppose that solid-body rotation was 
enforced in all convective regions. Another much more impor- 
tant is dynamical shear instabilities. These instabilities can en- 
sure that the rotation velocity is constant on equipotential sur- 
faces ( Pinsonneault et al. 1989) . 

The second category, with a time scale that is comparable 
to the Kelvin-Helmholtz time scale or the time scale for the 
evolution of the star, is that of secular instabilities. For secular 
instabilities, the transport process of angular momentum and 
chemical composition was treated as a diffusion process (Endal 
& Sofia 1978 ; Pinsonneault et al. ,1989) , so the radial equations 
for the redistribution of angular momentum and for the mass 
fraction X, are 



do. 



1 d 
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where is the diffusion coefficient. Because of some inherent 
uncertainties in the diffusion equation, the adjustable parame- 
ter fn is introduced to represent these uncertainties. Another 
adjustable parameter fc is used to account for how the insta- 
bilities mix material less efficiently than they transport angular 
momentum (Pinsonneault et al. ll989l) . The second term on the 
right-side of Eq. ST3[ is due to the secular contraction and/or 
expansion (Maeder & Zahn 1 19981 ), which can be dominant to 
induce differential rotation in a rotating star with a weak wind 
or without any wind. The second term on the right-side of Eq. 
(fT4b is the change in the nuclear reaction. The V, in the Eq. 
(fT4l) is the velocity of microscopic diffusion given by Thoul et 
al. (I1994 i In our model, we use the diffusion coefficient that 
was given by Zahn (I1993I I for secular shear instability. 



2c dlnT 



21G dr 



2dlnp^_^ I 
3 dr KpM(r) ' dr 



It has been thought that the Sun and T-Tauri with about 
one solar mass both possess an initial angular momentum of 
the order of 10^° g cm^s-^ (Kawaler ,1987. Paterno 1991) . 
The present total solar angular momentum is about lO'*^ g 
cm^.r '(Piipers |T998l Antia et al. |2000l Komm et al. "2007). 
Thus, the mean angular momentum loss rate is not less than 
10'*° g cOT^i^'yr"'. Using the formula for angular momentum 
loss (Mestel [T984l Kawaler [T988] l, 



dJ 2 dM , r^i , 
dt 3 dt ^ 



(16) 



and adopting the values of ^, and Q.^ as 10"**' g cm^s^'yr"', 
10"''* M^yr^^, and 10"^ rad s"', respectively, we can get the 
corotational radius r^ of solar wind is lOO/?©. It seems difficult 
for the solar wind to remain corotation with the Sun at a dis- 
tance of 100/?o- In this paper, we adopt a low initial angular 
momentum and ignore solar magnetic braking at the surface. 
This situation was also studied by Eggenberger et al. (I2005I I. 
Circulation is ignored because the circulation is very weak in a 
no-wind star, or it can even vanish in slow rotators (Zahn [T992T l. 
In our models, the total angular momentum is conserved, and 
the secular contraction and/or expansion are dominant to in- 
duce interior differential rotation. 



3.2. Magnetic effects 



Spruit ( 119991 ) reviewed the known magnetic instabilities in dif- 
ferentially rotating, stably stratified stellar interiors. Recently, 
Spruit (fT999l llOOa and Maeder & Meynet ([20031 12004] | de- 
veloped the Tayler-Spruit dynamo theory, which can generate 
magnetic fields in the radiative interior of differentially rotating 
stars. These fields are predominantly azimuthal components, 
B ~ Z?0. The existence of magnetic field and magnetic insta- 
bilities in the stellar interior provides a process in which angu- 
lar momentum and chemical compositions can be transported 
by magnetic stresses and magnetic instabilities. In this paper, 
we have assumed that the transport of angular momentum and 
mixing of material can be treated as a diffusion process. In the 
following, we give the diffusion coefficient for this process. 

In the plasma, the magnetic Prandtl number is 
(Brandenburg & Subramanian l2005l l 

)-'(— )-^ (17) 

^ ^ 20 



P,„ = - = 1.1 X 10-^(^)^( ^ 

7] 0.1 gcm- 

where v is the kinematic viscosity, rj the magnetic diffusivity, 
and the In A is the coulomb logarithm. In the solar radiative re- 
gion (r ~ 8 X 10'' - 2 X 10'' p ~ 20 - 0.2 gcm-^), the P„^ is 
about 10"^ - 10"^ The P,„ is about lO"'* - 10"^ in the solar con- 
vective zone (Brandenburg & Subramanian 2005). We assume 
In A - 10 in the above calculation. Thus, we ignore viscosity 
in the following equation. For a constant magnetic difiiisivity 
and shellular rotation, under axisymmetry and only consider- 
ing Lorentz force, the azimuthal components of the induction 
and momentum equations are 



dB^ ^ 1 
dt r^sirP-9 



(15) 11^ = 



do. 

'dt 



1 



Anpr^sirP-6 



, do. 
- = rsinOBr— , 

Bp ■ VirsinOB^) . 



(18) 



(19) 



If we assume that the effect of the magnetic diffusivity is to 
limit the growth of the toroidal field after some time (Mestel et 
al. I1987I I, the growth of the instability is halted by dissipative 
processes that operate on a timescale r, so that the second term 
on the left-hand side of Eq. ( fTSl l may be replaced simply by 
B^It (Barnes et al. 1999) . Substituting for the second term of 
Eq. ( fTST l and differentiating the Eq. ( fTST i with respect to time, 
we obtain (Barnes et al. 1 19991 ) 



d^B^ 1 dBi 



H — - rsinOB,- 

dh T dt 



do. 



(20) 
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For much longer times than the timescale of the instability, 
one would expect the term involving the first time derivative 
to dominate, so that (Barnes et al. 119991 ) 



dirsinOBt,) 



■V{rsin0B^)] 



(21) 



V ■ [D,nV{rsineB^)] . 



Equation ( 1211 1 is an approximate diffusion equation for rsinOB^ 
with a diffusion coefficient 



Dm — 



tB^ 



(22) 



For a steady equilibrium, the dissipating timescale has to 
match the growth timescale of the instability. In rotating stars 
under the condition cja « the growth time scale of the mag- 
netic instabihty is (Pitts & Tavler .1986.; Spruit , 1999i) 



Q 



where 



B 



lx>A 



is the Alfven frequency. If t- cr^^, one can get 



(23) 



(24) 



t^m — /!_„ 2 



(25) 



Note the diffusion coefficient D,,, derived from the induc- 
tion and momentum equations, so that one would expect that 
the angular momentum transport and the mixing of elements 
caused by magnetic fields obey a similar diffusion equation, 
such as Eqs. ( fT3] ) and (fT4l i with the diffusion coefficient D,,,. 
Although this diffusion coefficient is used in the equations 
of angular momentum transport and the mixing of elements, 
which is only an assumption, expression $25[ hints that this 
diffusion coefficient is related to the efficiency of the angular 
momentum transport. Expression ( |25] l can be rewritten as 



D,„ = r^Q^ . 



where 



Anp ' 

B^ 
Anp ' 



(26) 



(27) 



(28) 



4. Numerical calculation and results 

4.1. Choice of parameters 



The code originally written by Paczynski (|1969|[T970] l was up- 
dated by Sienkiewicz in 1995 and Yang et al. ( 120011 ). We mod- 
ified it to incorporate the hydrostatic effects of rotation on the 
equations of stellar structure, using the method of Kippenhahn- 
Meynet (Kippenhahn & Thomas 119701 Meynet & Maeder 
119971 ). The models are calculated using the OPAL equation of 
state (Rogers et al. '1996), OPAL opacity (Iglesias & Rogers 
119961 ), and the Alexander & Ferguson (1994) opacity table for 
low temperature. Element diffusion is incorporated for helium 
and metals (Thoul et al. 1994). The nuclear reaction rates have 
been updated according to Bahcall & Pinsonneault (119951 ). 
Energy transfer by convection is treated according to the stan- 
dard mixing length theory, and the boundaries of the convec- 
tion zones are determined by the Schwarzschild criterion. We 
adopt the solar age as 4.6x10** year, luminosity Lo=3. 844x10^^ 
ergi ', radius /?0=6. 96x10"^ cm, and the ratio of heavy ele- 
ments to hydrogen by mass Z/X=0.023 (Grevesse & Sauval 
fT998l ). 

As mentioned above, we assume that solid body rotation 
was enforced in the convective region. This assumption has 
been used by Pinsonneault et al ( 119891 ). Chaboyer et al ( 119951 ). 
and Huang (!2004bi The total initial angular momentum is a 
free parameter that is adjusted until the surface velocities of 
the solar-age models are near the solar surface velocity. We 
take the rotation rate of the Sun as a solid-body rotation, about 
2.72 xlO"^ rad/s (Komm et al. 120031 ). as the rotation rate of 
convective zone of our model. 

The strength and spatiotemporal distribution of magnetic 
fields inside the star are poorly known. Classical dynamo mod- 
els predict toroidal fields that are not stronger than about 10'* 
Gauss. But Choudhuri & Giknan (fT987] ). D'Silva & Choudhuri 
( 1993 ), and Caligari et al. ( 119951) have pointed out the value of 
the magnetic field at the bottom of convective zone as around 
10^ Gauss. The virial theorem (Parker' 1979) sets an upper limit 
to the magnitude of the average solar magnetic field: {B) < 10** 
Gauss. Dudorov et al. (119891 ) estimated the value of the poloidal 
magnetic field in the solar radiative zone to be about the order 
of unity. Using a kinematic model with prescribed internal ro- 
tation and a standard solar model. Fox and Bernstein (119871 ) 
investigated the existence of large-scale magnetic fields in the 
Sun, and found that the ratio of poloidal to toroidal components 
of the magnetic field is about 10 Although the ratio of B, to 
B may change with the radial coordinate r, the rotation rate Q, 
and the time f, the relation of % to r, Q., and t is unclear. As a 
first test, we take Br IB to be a constant, and take the order of 
the ratio of Br to B as 10"^ in the radiative region. 

In order to study the effects of magnetic fields, we construct 
three different types of solar models. These models are labelled 
as follows: 



is the Alfven velocity square. If the angular momentum is trans- 
ported by an Alfven wave, and the magnetic field B dominates 
as an azimuthal field, so that the larger the radio v^^/v^, the 
more efficient the angular momentum transport. 



- Ml: Model with no rotation and no magnetic field but with 
element diffusion; 

- M2: Model with rotation and no magnetic field but with 
element diffusion; 
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Table 1. Model parameters 



parameter 


Ml 


M2 


M3 


a 


1.797 


1.790 


1.790 


Yo 


0.2708 


0.2707 


0.2707 


iZ/X)o 


0.0276 


0.0277 


0.0277 


Y, 


0.2437 


0.2439 


0.2441 


(Z/X), 


0.023 


0.0231 


0.0232 




0.709 


0.710 


0.711 


fn 




1.0 


6.5x10-3 


fc 




0.02 


0.03 



Note.-a is the mixing-length parameter; io and (Z/X)q are the 
initial chemical compositions; Kj and (Z/X)s are the surface com- 
positions at an age of 4.6 x 10' years; Ri^^z is the radius at the base 
of the convective zone. 



- M3: Model with rotation and a magnetic field and element 
diffusion, but without including the effects of secular shear 
instability. 

Some parameters are summarized in Table [T] The mixing- 
length parameters a, Zq, Xq, /n, and fc are free parameters. 
The Xq, Zq and a are adjusted until the solar-age model has the 
values of the solar luminosity, radius, and {Z/X)s. Then Fq is 
determined by Y{) - I -Xq-Zq. The parameter /q was adjusted 
to fit the radial profiles of the angular velocity, which were got 
from helioseismology in the solar interior. Parameter /c was ad- 
justed to get the best radial profile of the sound speed. Finally, 
Ys, {Z/X)s and Rbc, are the results of calculation. All models are 
constructed by evolving a fully convective, pre-main-sequence, 
one solar mass model to the age of the present Sun. 

4.2. Transport of angular momentum 

We define the angular momentum density J(r, t) and total angu- 
lar momentum A(t) as 

7(r,f)= yp/Q(r,f), (29) 

A(f)= J(rj)dr, (30) 
Jo 

where t is time. 

In Figs. [T^ and[TJ), we present the rotation rates and their 
differentiation with respect to radius r as a function r at the age 
of 4.6 G years. The radial profile of the rotation rate of model 
M3 is flatter than of model M2. In our models, the radial gra- 
dient of angular velocity is negative. But the helioseismology 
has revealed that the gradient is positive at a low latitude and 
negative at a high latitude at the bottom of the solar convec- 
tive zone (Schou et al. 1998). At around r-QJ Rq in models 
M2 and M3, a sharp radial change in the rotation curve takes 
place. The thickness of the layer of sharp radial change in the 
rotation rate is about 0.036 Rq. Helioseismoloy also shows that 
that layer exists, and the thickness of the layer is around be- 
tween 0.02 Rq and 0.05 Rq (Christensen-Dalsgaard et al.[T99T] 
Koso vichev [19961 Basu & Anti a [T997al Basu ll997bl Co rbard 
et al. [T998l Corbard et al. [19991 Charbonneau et al. [T999l l. The 



rotation rate of model M3 is lower for r <0.4 Rq and higher for 
r >0.4 Rq than for model M2, which has the same total angu- 
lar momentum as model M3. When M2 and M3 have the same 
surface rotation rate, the rotation rate of M2 is higher than it is 
for M3 in whole radiative region. 

The diffusion coefficient of hydrodynamical model Dd is 
very small in the inner portion of the radiative region. It is 
ineffective to transport angular momentum outwards. But the 
diffusion coefficient of magnetic model D,„ is large enough to 
transport angular momentum from the core to the convective 
zone. That is to say magnetic fields are more efficient than the 
secular shear instabilities in transport angular momentum. The 
magnetic fields not only transport angular momentum outwards 
but strengthen the coupling between the radiative and convec- 
tive zone. Thus, under with the same total angular momentum, 
the rotation rate of model M3 is lower for r <0.4 Rq but is 
higher for r >0.4 Rq than the one for model M2 (FigH^); the 
gradient of the rotation rate of model M3 is significantly re- 
duced beneath the convective zone (Fig.[T]3). 

The angular momentum density J and its differentiation 
with respect to r of models M2 and M3, both with the same 
surface rotation rate at the age of 4.6 Gyr, are shown in Figs.[T]: 
and[TJl. The total angular momentum of model M3 is 1.91x1 0"^^ 
g cm^g-', which is consistent with helioseismic results (1.900 
± 0.0015)xlO'*'^ g cOT^i-i (Piipers il998j , 1.91x10'*^ g cm^g-^ 
(Antia et al. [2000l l. or (1.94 + 0.05)x lO'*'^ g cm^s-^ (Komm et 
al. 120031 . but that of model M2 is 2.1 x 10"^** g cm^^"', which is 
higher than helioseismic results. The angular momentum den- 
sity J shows a maximum 5.51 X 10" g cm .r' at r=0.35 Rq 
of model M3, and 6.1 x lO"'^ g cm .r' at r=0.34 Rq of model 
M2. Using the rotation rates obtained by inversions of helio- 
seismology, Komm et al. (I2003I I gave the maximum of J is 
5.5 X 10^^ g cm .r' for the Michelson Doppler Image (MDI) 
data, or 5.32x10^^ g cm s^^ for the Global Oscillation Network 
Group (GONG) data near r ^0.38 Rq. 

Model M3 is consistent with the Sun in the total solar an- 
gular momentum, the maximum of angular momentum density 
and the surface rotation rate. The magnetic fields are very im- 
portant in the angular momentum transport in radiative interior. 
In our model, the core is rotating faster than the rest of the ra- 
diative interior However, helioseismic results show the solar 
core is rotating more slowly than the rest of the radiative in- 
terior (Chaplin et al. '1999', Elsworth et al. '1995 ; Tomczyk et 
al. 1 19951 1, or else it is a solid-body rotation (Charbonneau et al. 
119981 1. Some mechanism of angular momentum transport may 
be missed, or the value of Br/B should be larger than the one 
used in the core in our model. A gravity wave may be the candi- 
date for angular momentum transport in the solar core, because 
the conditions for the known hydrodynamic instabilities occur- 
rence are not satisfied in a core that is rotating slowly (Spruit et 
al.QHS. 

4.3. Mixing and its effects on the solar structure 

The radial profiles of hydrogen mass fraction X are plotted 
in Fig. |2^. In Figs. |2}3, [2]:, and |2}l, we compare the density, 
adiabatic index, and squared sound speed of our models with 
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Fig. 1. Angular velocity Q., angular momentum density J, and their differentiation with respect to radius as a function of radius 
in the present-day Sun. (a) The solid line is model M3 with total initial angular momentum A - 1.91 x 10^^ g cm^s^^ . The 
long-dashed line and the dotted line show model M2 with A = 2.1 x 10^^ g cm^s^^ and A - 1.91 x 10"*^ g cm^s^^ , respectively, 
(b, c, and d) The solid line shows model M3 with A = 1.91 x 10"^^ g cm^i and the long-dashed Une is M2 with A = 2.1 x 10"*^ 

2 -1 

g cm s . 



the helioseismic inversion results given by Basu et al. (|1997cl 
I2OOOI 1. The mixing of magnetic instabilities causes more of an 
increase of the hydrogen abundance just beneath the convec- 
tive zone, but a decrease in the convective zone and near r=0.3 
Rq, than that of model Ml and M2. The increase of X between 
0.5 Rq and 0.7 Rq causes a decrease in the mean molecular 
weight /i and density p, so that the difference in density p in- 
creases in that region (Fig.|2]3). Since the squared sound speed 
oc r/yu, reducing yu must lead to an increase in c^. The differ- 
ence of squared sound speed, in the sense (c| - c^„oiigi}lc^„odei^ 
has a significant decrease, from 0.004 down to 0.0009 at r=0.65 
Rq (Fig.|2}l), whereas the effect is the opposite between 0.3 Rq 
and 0.5 Rq. In Fig.|2]l, there is a negative value, -1.06x10 of 
6c^ jc^ for models M2 and M3 at the bottom of the convective 
zone. That region is where the sharp radial change in X hap- 
pens due to the helium settling. The strong mixing leads to this 
negative value. The effect of mixing to the adiabatic index is 
very small, and the increase in X causes a small increase in the 
adiabatic index (Fig.|2j;). 



Although material mixing can account for the bump in the 
sound-speed's squared difference between 0.6 Rq and 0.7 Rq, 
the mass conservation causes a decrease in X, an increase in 
fi, and hence an increase in Sc^jc^ between 0.3 Rq and 0.5 
Rq. This change is not expected, but it must occur in a mix- 
ing model. 

5. Discussion and conclusions 

Parameter /n is significantly less than unity in our model. This 
could be a consequence of overestimating the ratio of B,. to B. 
The average ratio of the poloidal to toroidal field strength given 
by Fox & Bernstein (1 1987 J is 10"'°. If we take a lower value for 
Br IB, such as 10"*, instead of the previous one, we find that /n 
is 1 order of magnitude. For simplicity, /n and fc were taken as 
constant so the assumption might not be strictly correct. If we 
take the parameter as fnif, t), then, the value of fn(r, t) can be 
more than 6.5x10 in most regions and times. But the uncer- 
tainty of magnetic Br IB is larger than that of /n. We take /n as 
a constant. This is also a reason for low fa- The parameter fc is 




Fig. 2. (a) Radial distribution of hydrogen mass fraction, (b) Density difference between the Sun and the model, (c) Adiabatic 
index difference between the Sun and the model, (d) Squared sound-speed difference between the Sun and the model. All the 
differences are in the sense (Sun) -(Model). The solid line refers to model M3. The long-dashed line refers to model Ml. The 
dotted line refers to model M2 with total angular momentum 2. 1 x 10^^ g cin^g^K 



less than unity, too. In the magnetic model, the angular momen- 
tum can be transported by both magnetic stresses and magneto- 
hydrodynamical (MHD) instabilities, but only the MHD insta- 
bilities can mix material. So, the efficiency of material mixing 
must be less than that of angular momentum transport, i.e. low 

fc 

The magnetic field is more efficient than the secular shear 
instabilities in angular momentum transport outwards. Thus the 
rotation rate of model M3 is lower for r <0.4 R@ and higher for 
r >0.4 Rq than is model M2, which has the same total angular 
momentum as model M3. The model with the magnetic field is 
consistent with the heliosiesmic results in the solar total angu- 
lar momentum, the maximum of angular momentum density, 
and the surface rotation rate. 

In this paper, we obtained the diffusion coefficient of mag- 
netic angular momentum transport, which depends on the ra- 
dius, rotation rate, and magnetic fields. Using this diffusion co- 
efficient, we find that the angular momentum can be efficiently 
redistributed by magnetic fields. The radio of Br to B may be 
10"^ in the solar radiative interior. The coupling between the ra- 



diative region and convective zone can be strengthened by mag- 
netic fields, hence the sharp gradient of rotation rote (Fig.flJ)) is 
reduced at the bottom of the convective zone. In our model, the 
thickness of the layer of the sharp radial change in the rotation 
rate is about 0.036 Rq, and the rotation is almost uniform in the 
radiative region. There is a maximum, S.SlxlO"*^ g cin^g^^, of 
angular momentum density J at r=0.35 Rq for the model with 
a magnetic field. Material mixing associated with the angular 
momentum transfer leads to a change in the distribution of the 
chemical compositions, so that there is noticeable improvement 
in the profiles of the sound speed and density. 
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